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This paper describes a method of calculating stage parameters and 
flow distribution of axial turbines. The governing equations apply to 
space between the blade rows and are based on the assumption of 
rotationally symmetrical, compressible, adiabatic flow conditions. 
Results are presented for stage design and flow analysis calculations. 
Theoretical results from the calculation system are compared with 
experimental data from low pressure steam turbine tests. 


Modern steam turbines expand water vapor in 20 to 25 stages from a 
pressure of 3400 psia, 1050° F to about 0.80 psia, 94° F and handle thereby 
an 1800:1 increase in volume flow. Subdivision of the steam path into 
high pressure, intermediate pressure, and low pressure cylinders 2 makes 
the magnitude of pressure, temperature, and specific volume changes 
technically acceptable within each of the three units. The large' volume 
flows at the low pressure end arc handled by adding the appropriate 
number of low pressure turbines. 

We will discuss one stop in the steam path layout of such a low pressure 
turbine: the method of determining blade flow angles and calculating flow 
distribution in a sharply divergent, locally transonic flow field (ref. 1). 

Comparison of analytical predictions with low pressure test turbine 
results is the second subject of this paper. The calculation method is based 
on the well known “streamline curvature method” assuming rotationally 
symmetrical flow. The resulting circumferentially averaged information 
is used as a basis for the profile design. 


1 The author gratefully acknowledges the contributions of a great number of col- 
leagues and expresses appreciation to Mr. J. E. Fowler and to the General Electric 
Company for permitting this publication. 

2 A “cylinder” designates a separate housing. 
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DERIVATION OF WORKING EQUATIONS AND CALCULATION 
PROCEDURE 


The equations of motion and the continuity equation are transformed 
from a cylindrical coordinate system into a system with tin* stream 
function and the axial distance as independent variables. The meridional 
streamlines, along which the stream function is a constant, are approxi- 
mated by polynomials. These transformed equations and an assumed 
isentropic pressure density relationship in stream direction allow deter- 
mination of the axial velocity derivative and the pressure gradient in a 
radial direction at any axinormal plane (called “station”) outside a 
blade row. Derivation of the equations is given and one cycle of an 
iterative solution is discussed in this section. 

The Euler equations for steady state inviscous rotationally sym- 
metrical flow without body forces in a cylindrical coordinate system, 
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and the equation of continuity 
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are transformed into the \p,z plane. P designates the static pressure and 
p the fluid density. Coordinate system and velocity components are shown 
in figure 1. \p is a stream function defined by the differential expression 

d\p = 2irprwdr—2irpru'dz = \p r 'drp\p z -dz (3) 

Setting u = wa (a being the streamline slope in the meridional plane) and 
v = Y/r, 2ir = k, leads to the following set of equations in the \ p,z plane: 
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Partial derivatives are indicated by the coordinate subscripts. 
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Equations (4), (5), and (6) combined yield the Bernoulli equation in 
differential form along a streamline 
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Equation (4) is used for the calculation of the radial pressure dis- 
tribution by iteration. In order to do this, w : is determined from ('((nation 
(7) and equation (8), assuming isontropic expansion along a streamline 
and known streamline shape. 
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Equations (4) and (9), applying to the fluid motion at any axial 
station i, (fig. 2), are supplemented by the energy equation that is valid 
along any streamline (i/q = constant) and links up thermodynamically 
with station i— 1. 


Wi 2 =(Kl)W\^+(K2) 


2y, 

(T j— 1 1 ) Pi— 1 _ \Pi-J 


+ co, 2 (r i 2 -r 2 i _ 1 ) (10) 

lb is the total velocity relative to the blade row at which exit station ? is 
located. A'l and A’2 are the carryover and reaction coefficients, denoting 
fraction of upstream kinetic energy and reaction energy accounting for 
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Figure 2. — Grid for setup of calculation. 


the kinetic energy at station i, respectively, co is the wheel speed and 7 the 
specific heat ratio of the medium. 


The shape of the streamlines u and o z is 
polynomial curve fit such that 

3 approximated by piecewise 
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The coefficients of matrices A and C are functions of spacing of stations. 
Those of matrices B and 1) are functions of spacing and streamline slopes 
at upstream and downstream boundaries. It is the radius matrix. Only one 
matrix inversion is necessary at the beginning of iterations, oy is found by 
numerical differentiation. 

Equations (4), (9), (10), (11), and (12), equation of state or steam 
tables, and relations for velocity triangles are the basic equations used 
in the iteration at a station. This system can be solved either in a design 
mode or in an analysis mode for each station. 

The interchangability of modes permits steam path design to proceed 
through both new and available cascade components, besides the usual 
routine layout of all new stages or the analysis of an existing turbine at 
arbitrary flow conditions. 

Th(> sequence of calculations for a design iteration is sketched in the 
following: 

Assume station i is to be iterated. Having just passed station i— 1, gas 
conditions, velocities, and streamline shape are known as functions of 1 pj 
and 2 ,_i. Slopes ay,- and their derivatives ay,y, and oy.y are determined from 
respective matrices. Reaction and carryover coefficients may be calculated 
internally or given by input. 
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Next follows determination of w*,-,- (eq. (9)) pressure derivative 
(eq. (4) ) and the pressure distribution 



Pn is the known root pressure at station i. 

From the known pressure ratio and upstream data at station i— 1, 
total velocity W a relative to blade row at station i can be determined 
from equation (10), and the static enthalpy ha can also be calculated. 
Gas specific volume is either calculated from gas law or found from com- 
puterized steam tables. An alternate simplified calculation, including the 
wet region, is possible using gas laws and variable specific heat ratio by 
numerically differentiating steam table values. 

The axial velocity component, derived from local velocity triangles 
and given Ty/Tn distribution, is used to force 
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to approach a given value of r? ip — rj: oot by rotating the total velocity 
vector W a at the root such that the modified root axial velocity is given by 
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and by adjusting I\i accordingly. 

A new radius distribution, 



is finally calculated using a modified axial velocity, 
w i j moi = w ij + (wn moi -wa) DC 
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(DC = Damping Constant.) 

Convergence is checked by comparing radii, slopes <r, and slope deriva- 
tive based on the new radii with those computed in the previous 
iteration. Depending on these checks, the described iteration is either 
repeated or the system moves on to the next station. This process goes on 
until the downstream boundary is reached. The system will either return 
to the upstream boundary for a new iteration loop or move to answer 
calculations, depending on an overall streamline radius convergence 
criterion. 
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NUMERICAL EXAMPLES 
Design Calculation 

In a design case, the whole low pressure turbine is calculated in one or 
two sections. In the latter case an overlay is made 1 at the joining stations 
in order to ensure streamline as well as thermodynamic continuity. Design 
information for one stage out of several calculated is shown in figures 3 
through 6. 

Figure 3 shows sin a (a is the nozzle exit angle) versus radial height 
and the partial derivative of sin a with respect to z at constant radius. 
The latter information allows one to account for nonradial trailing edges 
and other small axial adjustments during the course of design. Figure 4 
relates the relative bucket entrance angle and its derivative with respect 
to z, and figure 5 describes the sine of the relative bucket exit angle y and 
its derivative with respect to z as a function of the radius. Finally, in 
figure 6, the absolute stage exit angle versus radius is given. This in- 
formation is the basis for the next step in designing the blade sections. 
Additional information such as pressure, velocity, .Mach number dis- 
tribution, and streamline shape is available, serving as a guide in the 
successive progress of a low pressure turbine layout. 

Analysis of Test Conditions and Comparison of Theoretical Data 
with Laboratory Measurements 

Calibration of the outlined procedure is of most interest. The cal- 
culation system is used for this purpose in its analysis mode to make 
predictions of pressure distributions, flow distributions, and flow angles 
at any station based on given turbine geometry, estimated efficiency, root 



Figure 3. — Design information: sine of nozzle exit angle and its derivative versus radius. 
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Figure 4. — Design information: relative bucket entrance angle and its derivative versus 

radius. 



Figure 5. — Design information: sine of relative bucket exit angle and its derivative 

versus radius. 

stage exit pressures, upstream and downstream conditions, and total flow. 
These predictions are compared with available traverse data behind the 
second, third, and fourth stage of a test turbine. Since most of the test 
data were accumulated under off-design conditions due to high initial 
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Figure 6. — Design information: absolute stage leaving angle. 


superheat necessary for temperature traverses, the comparison makes for 
an especially good test of the capabilities of the analysis mode. 

The test data presented here was obtained in the low pressure turbine 
test facility in the Product Development Laboratory of the Large Steam 
Turbine Department of General Electric Company. A detailed description 
of this facility is given in reference 2. The last four stages of a 30-in. low 
pressure turbine were tested. A cross section of the test turbine showing 
location of pressure taps (P), and temperature sensors (T), as well as 
traverse locations behind the second, third, and fourth stage is exhibited 
in figure 7. Figure 8 is a photo of the test rotor. 

Test data (circles) and theoretical predictions (full lines) are plotted 
versus the respective radius ratio in figures 9 through 13. 

Figures 9 and 10 present traverse data behind exit from stage number 
two. There is maximum discrepancy in the static pressure readings of 
0.18 psia or a 2.5 percent deviation in absolute pressure level. The absolute 
leaving angle is negative due to lower than design velocity ratios. Cal- 
culated angles are too large, on the average by 7.8°, when the effect of tip 
leakage 1 is discarded. The deviation in both graphs is in the same direction. 

Comparisons of conditions at the third stage exit are shown in figures 
11 and 12. The measured axial velocity decreases towards the tip less 
rapidly than indicated by the analysis. The leaving angle is again nega- 
tive, indicating that the direction of the tangential leaving velocity is 
opposite to direction of wheel rotation. Calculations show angles larger 
on the average by 2.5°. 

Average carryover and reaction coefficients were used in the cal- 
culations of the second and third stage results. In the analysis of the last 
stage, an estimated nozzle efficiency and the known stage streamline 
efficiency were used for an approximation of the bucket efficiency. The 
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Figure 7. — Cross section through test turbine. 


graph in figure 13 compares predicted and measured axial velocity dis- 
tribution on a dimensionless basis. 

Finally, we want to call attention to the different flow regimes that do 
occur in different portions of a large radius ratio low pressure steam 
turbine stage. The absolute velocity entering the stator is well subsonic, 
while the stator exit velocity ranges from supersonic at the root to sub- 
sonic at the tip as shown by the Mach number plot in figure 14. The 
relative exit velocities at the rotor change from transonic velocities near 
the root to supersonic velocities along the outer portion of the bucket. 
The calculation system discussed is able to analyze these stages. 
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Figure 8. — Thirty-inch LSB test turbine rotor. 
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Figure 9. — Comparison of predicted and 
measured static pressure distribution at exit 
of stage 2. 


igure 10. — Comparison of predicted and 
measured absolute leaving angle at the exit 
of stage 2. 




Figure 11. — Comparison of pre- 
dicted and measured axial ve- 
locity distribution at the exit of 
stage 3. 
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Figure 12. — Comparison of predicted and meas- 
ured absolute leaving angle at exit of stage 3. 



ABSOLUTE LEAVING ANGLE IN DEG. 
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Figure 14. — Mach number dis- 
tribution at last stage nozzle 
exit. 


CONCLUSIONS 

The system described has proved to be a reliable design tool. Its simple 
structure allows updating of the system by input of experimental data 
derived from tests in air and steam. The analysis mode is helpful in 
interpretation of test evidence. Comparison of analysis with test data 
shows that the system overshoots traverse data at the second stage 
traverse, while third and fourth stage traverses are predicted well. The 
match of static pressure and leaving angle depend very much on the 
effective nozzle areas. Differences in the axial velocity distributions can 
be reduced by taking account of radial bucket efficiency distribution. 
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DISCUSSION 


G. K. SEROVY (Iowa State University) : The author has presented 
a method for solution of a most difficult problem. The basic approach, 
as indicated by the author, is not new. It has been used to attack both 
design and analysis problems in axial-flow compressors, as reported inde- 
pendently by Smith and by Novak in ASME papers several years ago. 
Another closely related method was described by Renaudin and Somm 
of the Brown Boveri Company at a Symposium on Flow Research on 
Blading held in Baden, Switzerland, in 19G9. The authors of earlier 
papers as well as this one have all succeeded in handling a very compli- 
cated situation in numerical analysis which has been encountered by all 
those who attempt to iterate on stream-surface shapes, especially when 
the flow is of a transonic character. 

It would be helpful in evaluating the paper and results if the author 
could tell us something about the initial conditions assumed at the up- 
stream boundary. Also he must have in his solution system some method 
for making an initial guess at the shapes and locations of the (\p = constant ) 
lines. How is this done? 

In the case of axial-flow compressor calculations, it has been possible 
to improve solutions for the design case by using stations located inside 
the blade rows. It has always seemed probable to this observer that the 
blade-to-blade flow problem in turbines has been under control, so that 
such solutions (intra-blade stations) might be very feasible for steam 
turbines. 

One might also take a certain amount of fiendish pleasure in noting 
that the “turbine crowd” has a most serious problem in accounting for 
end wall flows. We have this problem in equal or greater quantity in the 
“compressor crowd” and it is to be hoped that we can cooperate in de- 
veloping consistent procedures for doing a better job in these regions. 
Mr. Herzog deserves our thanks for giving us a clear outline of his method, 
supported by comparison with experimental radial surveys. 

H. D. LINHARDT (Airco Cryogenics) : I understand the turbine 
you have discussed is operating under wet steam conditions in the last 
stage. Could you please define what are the wet steam conditions in the 
last stage, what is the percentage of wetness and how does the performance 
change when operating under wet conditions? It would be interesting to 
know what kind of design procedure you use and what the condensate 
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droplet sizes are in the last stage. How do you correlate performance and 
erosion phenomena with droplet size? 

R. M. HEARSEY (Ohio State University Research Foundation) : I 
wonder if you have considered modifying your procedure slightly from 
what I understand it to be to include the possibility of the calculation 
stations being nonradial, which would appear to allow you to calculate 
velocity distributions rather closer to the blade edges. 

C. EARN (Westinghouse Research and Development) : What kind 
of loss criteria did you use in the computation, especially for the last 
stage? 

You only use two axial stations per stage. I always wonder whether 
this is accurate enough. 

Is this a design program or performance program? If it’s a design 
program, what design parameters are you specifying? 

HERZOG (author) : The initial conditions assumed to be given at the 
upstream boundary are the axial location of the station, the stagnation 
pressure and stagnation enthalpy, as well as the radius and slope at each 
streamline. Total conditions can vary from streamline to streamline. 

The initial guess at shape and location of the (^ = constant) lines is 
made by subdividing the flow path into a given number of equal area 
filaments at each station. This provides a complete radius matrix that is 
used to determine initial matrices for streamline slope <r,y and slope 
derivative <r zi j. 

Introduction of intra-blade stations is desirable and expected to im- 
prove accuracy of shape and location of the {\p = constant) lines. As of 
now, we use only three stations per stage. Additional stations can be 
introduced into the vaneless portions of a stage. Use of nonradial stations 
is not planned at this time. Information needed along nonradial contours 
such as blade trailing edges are determined from information at the 
radial station which is modified by linear extrapolation either in streamline 
or in axial direction. 

The program is mainly used as a design tool for large radius ratio 
steam turbines; we specify, in addition to upstream conditions, the inner 
and outer wall shape, blade spacing in axial direction, axial energy and 
reaction distribution at blade roots, and the stage total flow, wheel 
speed, and pressure distribution at the turbine exit. Existing components 
may be incorporated in new designs. 

Steam entering the last stage under discussion is about 3§ percent 
wet and leaves the stage at about 10 percent wet. 

The performance of the turbine deteriorates with increase in initial 
moisture content. This effect is shown in figure D-l, which was taken 
from reference 2. This graph displays the change in turbine efficiency 
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with the change of the initial moisture content at the inlet of a four 
stage low pressure tost turbine. In the range of super-heated inlet con- 
dition only part of the turbine is running on wet steam. We did not 
measure drop sizes in the turbine. Size measurements in a shock tunnel 
indicated an initial drop dimension, after spontaneous condensation took 
place, of submicron size. It is believed that these small drops do not 
cause erosion and that most of the damage is done by much larger “sec- 
ondary drops” torn off the partition trailing edges. We are not aware of 
any laws correlating performance and erosion phenomena with drop size. 




